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In this paper, a time variable cosmological constant (CC) from renormalization group equations 
(RGEs) is explored, where the renormalization scale /j 2 = R~qq = Max(H + 2H 2 , —H) is taken. The 
cosmological parameters, such as dimensionless energy density, deceleration parameter and effective 
equation of state of CC etc, are derived. Also, the cosmic observational constraints are implemented 
to test the model's consistence. The results show that it is compatible with cosmic data. So, it 
' would be a viable dark energy model. 
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I. INTRODUCTION 



CO 

Cosmological constant (CC) is a long standing issue in cosmology and physics. It was first introduced by Einstein 

u: 

to realize a static universe about a century ago. However, it was found that this space-time was unstable and not 
consistent with observed cosmological expansion. Recently, CC returns to cosmology as a natural candidate to dark 

O . energy to explain recent cosmic observation that our universe is undergoing an accelerated expansion firstly deduced 
from observational results of type la supernovae 0,0. In the context of quantum field theory (QFT), CC has relations 

r*, \ with the vacuum or zero point energy density of quantum fields, via 

£3 . 

£ 1 1 f A ^k 2 dk fTo o A4 

where A 3> m is a UV cut-off. To balance an assumed UV cut-off A and the observational smallness of CC, tremendous 
fine tuning is required. This is the so-called cosmological constant problem Q. As known, renormalization in QFT 
can handle the infinities, and a dependence of the renormalized constants on some energy scale /z is leaded. This 
renormalized scale is usually identified with external momentum or characteristic scale of the environment. 

QFT in curved space-time leads to an infinite effective action or vacuum expectation values (VEV) of the energy- 
momentum tensors of the fields. A renormalization treatment can yield a scale-dependent or running CC and a running 
Newton constant. Though, the absolute values can not be calculated, the change with respect to the renormalization 
scale can be calculated via RGEs originating in QFT Q and quantum gravity ||. In the cosmological context, it is 
reasonable to identify the renormalization scale with some characteristic scales of cosmology. The related investigation 
can be found in where the renormalization scale /i was given by the Hubble scale H, the inverse radius R^ 1 of 
. £^ ' the cosmological event horizon and the inverse radius T _1 of the particle horizon. A scaling law having decoupling 
behavior at low energy in H (which was extensively studied in 0]) was explored where the corresponding /3- function 
$— i for A is 



H— = V, (2) 
a/j, 

where A ~ ±M 2 is determined by the masses M and the spins of the fields. Assuming constant masses and [iq ~ Hq, 
one has @ 

A(/x) _ [l 2 M 2 

— — -L + Li — , Li ~ ±— j, (3) 

A Mf, 

where Lo = 1 — L\. For sub-Planckian masses M, \Li\ -C 1. Then Lq ~ 1 as a model parameter can be teste by 
cosmic observations. Obviously, when Lq = 1, the CC becomes scale independent and a real CC is recovered. 
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Another consideration of time variable CC can be found in @, H, H, OH [Tl[. Following the work @, we are 
going to take Rq 2 c — Max(H + 2H 2 ,—H) as a possible candidate to a renormalization scale to investigate the 
evolution of our universe in this work when the Newton constant is fixed. Though, the reports of [13] have shown 
that the cosmology is not a RG flow where RG was checked by a massless, minimally coupled scalar with a quartic 
selfinteraction on a nondynamical, locally de Sitter background. But it does not preclude using the RG conventionally 
to relate quantities at different constant scales. So, R^c can potentially be used as a renormalization scale. Here, 
R^c = Max(H + 2H 2 , —H) is the causal connection scale for spatially flat universe [HI]. That was investigated in 



;iaiiy ne 

the context of holographic dark energy [13l | and Ricci dark energy 2 Furthermore, it was found that only the 
case where Rq 2 c = H + 2H 2 as an IR cut-off was consistent with the current cosmological observations when the 
vacuum density appears as an independently conserved energy component [HI . However, when a time variable CC is 
considered these two cases must be checked over again for its coupling with cold dark matter. It is just the case that 
will be explored in this paper. 

This paper is structured as follows. In Section HH we give a brief review of time variable CC. We will discuss the 
evolution of time variable CC and derive the cosmological parameters in Section Hill Cosmic observational constraints 
are implemented in Section HVl Section [V] are conclusions. 

II. TIME VARIABLE CC 



The Einstein equation with a cosmological constant is written as 

1 
2 



R^u - t: r 9^ + A 6V = SivGT^, (4) 



where T M „ is the energy-momentum tensor of ordinary matter and radiation. From the Bianchi identity, one has the 
conservation of the energy-momentum tensor V^Ij,,, = 0, it follows necessarily that A is a constant. To have a time 
variable cosmological constant A = A(i), one can move the cosmological constant to the right hand side of Eq. (|4|) 
and take = T M „ — ^^g^u as the total energy-momentum tensor. Once again to preserve the Bianchi identity or 
local energy-momentum conservation law, V^T^j, = 0, one has, in a spatially flat FRW universe, 

Pa + Pm + 3i? (1 + w m ) p m = 0, (5) 

where p\ = MjA(i) is the energy density of time variable cosmological constant and its equation of state is wa = —1, 

and w m is the equation of state of ordinary matter, for dark matter w m — 0. It is natural to consider interactions 

between variable cosmological constant and dark matter as seen from Eq. After introducing an interaction 
term Q, one has 

p m + 3H (1 + w m ) p m = Q, (6) 

PA + 3H (pa+Pa) = -Q, (7) 
and the total energy-momentum conservation equation 

ptot + 3H ( Ptot + p tot ) = 0. (8) 

For a time variable cosmological constant, the equality pA + Pa = still holds. Immediately, one has the interaction 
term Q = —Pa which is different from the interactions between dark matter and dark energy considered in the 
literatures [16| where a general interacting form Q = 3b 2 H (p m + pa) is put by hand. With observation to Eq. (JT)), 
the interaction term Q can be moved to the left hand side of the equation, and one has the effective pressure of the 
time variable cosmological constant- dark energy 

PA + iH[p A +PA ff ) =0, (9) 

where p e j[* = pa + Mj is the effective dark energy pressure. Also, one can define the effective equation of state of 
dark energy, (for other definition, please see [13), 



w 



A 



■■ff _ Ea_ 

PA 

-1 + 



Q 



3Hp A 



1 d In pa . 
3 dm a 
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The Friedmann equation as usual can be written as, in a spacially flat FRW universe, 

H 2 = TTW (Pm + PA) ■ (11) 



III. EVOLUTION OF TIME VARIABLE CC AND COSMOLOGICAL PARAMETERS 

The time variable CC and G was explored in [g], where the renormalization scale fi was given by the Hubble scale 
H, the inverse radius Rr 1 of the cosmological event horizon and the inverse radius T _1 of the particle horizon. In this 
paper, we are going to reconsider the time variable CC when R~^ c = Max(H + 2H 2 , —H) is taken as a renormalization 
scale say p? = Rqq- Though the authors of [l3| have claimed the case where Rqq = H + 2H 2 as an IR cut-off was 
consistent with the current cosmological observations when the vacuum density appears as an independently conserved 
energy component. For the existence of effective interaction between time variable CC and cold dark matter, the two 
cases must been checked over again. We name the case p, 2 = H + 2H 2 Model A and p 2 = —H Model B. 



A. Molde A: p? = H + 2H 2 

In this case, the time variable CC can be written as 

R 

A(f) = A (L + ii^) 

_ H + 2H 2 . 



(12) 



From the above equation, it seems the CC behaves quite like Ricci dark energy. But in fact, it is different from that 
for its effective interaction with cold dark matter. The corresponding vacuum energy density is 



Pa = M 2 A(t) 



Af£Ao(£o + £i 



H + 2H 2 



= Mj,A 



H + 2H 2 ' 
L Q + M Q {H + 2H 2 ) 



(13) 



where Mo = (1 — Lq)/(Hq + 2Hq). For its interaction between p\ and cold dark matter p m , they are not con- 
servative separately. By using the definition of dimensionless density parameters ft m (z) = p m /(3MpH 2 ) and 
tt A (z) = p A /(3M],H 2 ), one has 



i-n A 



PA- 



Then after easy algebra, the conservation Eq. ((5]) can be rewritten as 

h + \h 2 {\ - n A ) = o. 



(14) 



(15) 



From Eq. (JT3J), one has the expression of 17/ 



n A = 



Ao 
3H 2 



L + M (H + 2H 2 ) 



Ao 2 . 
3H 2 3 



H 
H 2 



(16) 



where Aq — AqMq and Bq = AqLq for convenience. Combining Eq. (|15p and Eq. one has 



' H 2 2 //-' ' 2 1 



2Ac 



= 0, 



(17) 
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In terms of redshift z, the above Eq. (fTTf can be rewritten as 

A \ (1 + z) d\nE(z) B' Q 1 



dz 



3 / 2Ap 
2 E(z) ' 2 I 3 



0. 



(18) 



where H 2 (z) = HqE(z) is used and B' = Bq/Hq which is a dimensionless parameter. The above differential equation 
has the integration 



E(z) = 



1 - 



B' 



(3 - 2A ) 



(1 + z f+Ao/{A -2) + 



B' 



(3 - 2A ) ■ 



(19) 



B. Model B: fi 2 = —H 

In this case, the time variable CC can be written as 

A(t)=A (L + L 17 -). 

Ho 



As done in IIII A| one has the expression of Q\ and differential equation of H 



n A = 



Ao 
3H 2 



L + N H 



[ 2 'h 2 2 H 2 2 ' 

where N = (1 — L )/H and Co = A A . Also, one can find the Hubble parameter H(z) as a solution of Eq. 
with respect to redshift z as follows 

H 2 (z) = H 2 {(l-^) (l + 2 )^ + |« 
It is clear that ACDM is recovered when Cq — 0, i.e. Lq = 1. 



(20) 

(21) 
(22) 



(23) 



Discussion 



From these Friedmann equations ([19]) and Eq. (|23|) , one can immediately find out that the first terms of the right 
hand of the equations behave like cold dark matter for A = and Co = respectively, i.e. L = 1. In these cases, 
the ACDM universe are recovered as expected in introduction. These models contain two parameters ^4o(Co) an d 
B' Q which can be determined by cosmic observations. If this model does not badly depart from ACDM universe, we 
can estimate the values of parameters A (Co) and B' Q . It is to say A ~ 0(Co ~ 0) and B' ~ 3Qao ~ 2.1 which can 
be tested by cosmic observations. In terms of redshift, the deceleration parameter and effective EoS of CC can be 
written as 



-1 



w 



eff 



{l + z)d\nH 2 

2 dz~ ' 
(l + z) fd\nH 2 d1nSl A (z) 



dz 



dz 



(24) 
(25) 



From Eq. (|19p . one can find two singularity points with parameter values of A$ = 2 and 3/2 respectively. The same 
case can be found in Eq. (f2"3"]) when Co = 2. When Ao = 2, the first term of left hand side of Eq. (TTT]) vanishes. Then 
one has a constant Hubble parameter H 2 = —1/Bq. For the positivity of the value of B > 0, it does not describe a 
physical system. However, for Co = 2, one has H 2 = -Bo/3. Then, in this case, a de Sitter or anti de Sitter universe 
can be obtained. When A — 3/2, one has H = 2B which corresponds to scale factor a(t) ~ exp(_B ^ 2 )- 
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Model 


Xmin 


^Ao(lcr) 


A or Go (Iff) 


Bo(lff) 


Zr(lff) 


XLin/dof 


A 


473.593 


0.732 ± 0.021 


(1.0 x 10- 11 )™ 7 


2 1Q6+ 0780 
z,.iitu_ 0833 


762+ 071 
u. ( uz_ ogg 


1.193 


B 


473.593 


0.732 ±0.021 


(-7.011 x lO- 10 ^™ 


2 1Q6+ 0780 

*- Lao _0 .0833 


762+ ' 072 


1.193 



TABLE I: The minimum values of \ 2 an d best fit values of the parameters. Here zt is the transition redshift from decelerated 
expansion to accelerated expansion and do} denotes the model degrees of freedom. 



IV. COSMIC OBSERVATIONAL CONSTRAINTS 

Now, it is proper to present the constraint results by using cosmic observations: SN la, BAO and CMB shift parame- 
ter R, for the details please see Appendix [Al In this work, 397 SN la Constitution dataset, the ratio Dy(0.35)/Z)y(0.2) 
detected by BAO and CMB R from WMAP5 are used. After the calculation as described in Appendix fAl the results 
are listed in Tab. Q] The evolution curves of q(z), to^ ( z ) anc ^ dimensionless density parameters Vt m [z) and £l\(z) 
are plotted in Fig. Q] and Fig. [2] Also, the contour plots of Aq — B' and Co — B' are shown in Fig. [3] and Fig. [4] 




FIG. 1: Model A: The evolution curves of q(z) (left panel), w^' (z) (central panel) with 1<t error region and dimensionless 
parameters fi m (z) and Qa(z) (right panel) with respect to redshift z where the best fit values are adopted. 




FIG. 2: Model B: The evolution curves of q(z) (left panel), w A (z) (central panel) with la error region and dimensionless 
parameters Q m (z) and Qa(z) (right panel) with respect to redshift z where the best fit values are adopted. 

From the left panels of Fig. [1] and Fig. [21 one can easily find that our universe is undergoing accelerated expansion at 
late time, and the transition redshift from decelerated expansion to accelerated expansion are zt ~ 0.762, which are 
consistent with other analysis results with best fit parameter values. In the early epoch, our universe is dominated 
by cold dark matter, that can be seen from the left and right panels of Fig. [1] and Fig. [5] where the deceleration 
parameter is q — > 1/2 (dark matter dominated) and f2 m — * 1 at high redshift. From the central panels of Fig. [1] 
and Fig. [2l one can see the effective EoS of time variable CC is almost constant te^ \z) ~ — 1. So, the universe is 
quasi-ACDM and it is not necessary to worry about the structure formation of the universe. With the best fit values 
of the parameters and Lq = B' / (3J7ao), one obtains the mass of the fields is M = |1 — L^^Mp ~ 1.490 x 10 _8 Mp 
and L\ ~ —2.220 x 10 -16 in Model A and Model B which confirms the prediction in the introduction. 

V. CONCLUSIONS 

In this paper, time variable cosmological constant from renormalization group equations (RGEs) is explored, where 
the renormalization scale /i 2 = R^c — Max(H + 2H 2 ,—H) is taken. The cosmological parameters, such as dimen- 
sionless energy density deceleration parameter and effective EoS of CC etc, are derived. Also, the comic observational 
constraints are implemented to test the model's consistence, the results are shown in Tab. [U As investigated, with 



FIG. 3: Model A: The contours in the planes of Aq — B' with la and 2a regions. The central dot denotes the best fit values 
of model parameters. 



FIG. 4: Model B: The contours in the planes of Co — B' with la and 2a regions. The central dot denotes the best fit values 
of model parameters. 

this time variable CC, the universe is undergoing an accelerated expansion at late time and a decelerated expansion 
at high redshift. And, the transition redshift from decelerated expansion to accelerated expansion is zt ~ 0.762 which 
is consistent with other results. The effective EoS of time variable CC is almost constant w e J^{z) ~ —1. In early 
epoch, our universe is dominated by cold dark matter, that can be seen from the left and right panels of Fig. [T] and 
Fig. [5] where the deceleration parameter is q — > 1/2 and H m — ► 1 (dark matter dominated) at high redshift. And the 
current cold dark matter density ratio is Oao ~ 0.732 which is also compatible with other analysis. So, via RGEs 
with the renormalization scale /j 2 = Rcc = Max(H + 2H 2 , —H), time variable CC is viable. With the best fit values 
of the parameters and Lq = Sq/(30ao)j one obtains the mass of the fields is M = |1 — Lq^'^Mp ~ 1.490 x 10 _8 Mp 
and L\ ~ —2.220 x 10~ 16 in Model A and Model B which confirms the prediction in the introduction and implies 
that the mass scale having relations with CC about 10 -8 order of Mp. 
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APPENDIX A: COSMIC OBSERVATIONS 



1. SN la 



We constrain the parameters with the 397 SN la Constitution dataset including 397 SN la [18j . Constraints from 
SN la can be obtained by fitting the distance modulus 



Hth(z) = 5log 10 (D L (z)) + /i , 
where, Dl(z) is the Hubble free luminosity distance Hq<1l(z)/c and 

dz' 

d L (z) = c(l + z) ' 



(Al) 



/o ™ 
Ai = 42.38-51og 10 /i, (A3) 

where H is the Hubble constant which is denoted in a re-normalized quantity h defined as H — 100ft. km s _1 Mpc _1 . 
The observed distance moduli Hobs(zi) of SN la at Zi is 

Vobs(zi) = m obs (zi) - M, (A4) 

where M is their absolute magnitudes. 

For SN la dataset, the best fit values of parameters in a model can be determined by the likelihood analysis is based 
on the calculation of 

2/ s IVobs(Zi) - HthiPs.Zi)] 2 



SNIa 



07 



E 

SNIa 



[5 log 10 (D L (p s ,Zi)) - m obs (zi) + mo] 



(A5) 



where mo = /Jo + M is a nuisance parameter (containing the absolute magnitude and Hq) that we analytically 
marginalize over [l9j . 



X 2 (p s ) = -21n 



to obtain 



exp 



B 



1 



r = A - — + In — 



X 2 (Ps,m ) 



C 



dmo 



(A6) 



C 



2tt 



where 



A = E 



[51og 10 (Di(p 8 ,^)) -m bs(zi)Y 



SNIa 



^=E 



5log w (D L (p s , Zi) - m obs (zi) 



SNIa 



(A7) 

(A8) 
(A9) 
(A10) 

(All) 

is used instead of Eq. (|A7|) to perform the likelihood analysis. They arc equivalent, when the prior for tuq is flat, as 
is implied in (|A6[) . and the errors Ui are model independent, what also is the case here. 

To determine the best fit parameters for each model, we minimize % 2 (p s , B/C) which is equivalent to maximizing 
the likelihood 



SNIa Qi 



The Eq. (|A5[) has a minimum at the nuisance parameter value rrio = B/C. Sometimes, the expression 

X 2 S Nia(Ps,B/C)=A-(B 2 /C) 



C( Ps ) OC e -X 2 (p s ,B/C)/2^ 



(A12) 
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BAO 



The BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, and measure 
the distance-redshift relation at z = 0.2. BAO in the clustering of the SDSS luminous red galaxies measure the 
distance-redshift relation at z = 0.35. The observed scale of the BAO calculated from these samples and from the 
combined sample are jointly analyzed using estimates of the correlated errors, to constrain the form of the distance 
measure D v (z) [13, HE H3] 



D v (z) 



{l + zfD A {z)- 



1/3 

(A13) 



H(z) 

where Da(z) is the proper (not comoving) angular diameter distance which has the following relation with (Ll(z) 

Matching the BAO to have the same measured scale at all redshifts then gives [22| 

Dv(0.35)/Dv(0.2) = 1.812 ± 0.060. (A15) 

Then, the X%ao(Ps) is given as 

x1aM _ im/D^)-^ (A16) 

3. CMB shift Parameter R 

The CMB shift parameter R is given by [23[ 

R(z*) = yJn m H§{l + z*)D A (z*)/c (A17) 



which is related to the second distance ratio Da (z*)H (z*)/c by a factor y/1 + z*. Here the redshift z* (the decoupling 
epoch of photons) is obtained by using the fitting function [24j | 



Z-k 



1048 [1 + 0.00124(fi 6 /i 2 )- 738 ] [1 + gi(n m h 2 ) 92 ] , (A18) 
where the functions gi and g2 are given as 

,9i = 0.0783(^/i 2 )~ - 238 (l + 39.5(r! b / l 2 ) - 763 ) _1 , (A19) 

g 2 = 0.560 (1 + 21. l(n 5 /i 2 ) x ' 81 ) -1 . (A20) 



The 5-year WMAP data of R(z*) = 1.710 ± 0.019 [2J| will be used as constraint from CMB, then the Xcmb(Ps) is 
given as 

2 (R(z*) - 1-710) 2 
Xomb(P») = FT7j77j2 ■ ( A21 ) 
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